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1. Introduction 


Quantum Field Theory (QFT) models are usually defined by specifying the La- 
grangian of the theory. Such definition is not unique, in the sense that it allows for 
transformations of the field basis, leading to equivalent descriptions of the model, with 
different degrees of freedom. The calculation of transition amplitudes can be performed 
in any basis, however two special choices are most convenient from the practical point 
of view. 

In many cases the Lagrangian of the model is initially constructed in terms of fields 
having definite charges under some symmetry groups - we call it symmetry or interaction 
basis. However, in general, fields defined in such a way do not correspond to physical 
degrees of freedom of the theory, and some redefinitions have to be performed in order to 
identify the physical fields (with the spontaneous symmetry breaking being the typical 
example). 

Another possibility is to use mass eigenbasis, in which bilinear terms (kinetic and 
mass matrices) in the Lagrangian have been diagonalized and states of the theory corre¬ 
spond to physical particles with definite mass. The transformation from the initial basis 
to the mass eigenstates basis is performed by unitary rotations (“mixing matrices”) in 
the field space. Perturbative calculations of the amplitudes in the mass eigenstates ba¬ 
sis lead to results expressed in terms of physically measurable quantities, i.e. physical 
masses and the elements of the mixing matrices. They are usually more compact as 
compared to those obtained in any other basis and best suited for numerical computa¬ 
tions. However, the analytical dependence of such amplitudes on the initial interaction 
basis is typically complicated and difficult to use for qualitative interpretation. 

For the latter purpose, it is often useful to have analytic expressions for the tran¬ 
sition amplitudes calculated directly in the interaction basis. They can be obtained 
using two different methods. Firstly, one can perform an independent diagrammatic 
calculation of the amplitude using the aOTroximation commonly referred to as the Mass 
Insertion Approximation (MIA) i, 3 3. In this approach, diagonal elements of the 
mass matrices are absorbed into the definition of (unphysical) massive propagators and 
the amplitude is, at every loop order, expanded in an infinite series of non-diagonal el¬ 
ements of mass matrices, commonly referred to as mass insertions (MI). Alternatively, 
as proven in Ref. 3 > fhe MIA result can be obtained directly from the mass eigenstates 
amplitude employing the purely algebraic technique coined in Ref 3 the “Flavor 
Expansion Theorem” (FET)^. The last method, the FET expansion, has two important 
advantages. Firstly, it allows to avoid the Feynman diagram calculation with mass in¬ 
sertions, which is usually tedious and prone to errors or omissions of important terms. 
Secondly, it can become, to large extent, automatized. 

In this paper we describe the MassToMI package, written with the use of Mathemat- 


^The first non-trivial order of the FET expansion, with applications in MSSM flavor physics, have 
been presented in refs. Higher orders could be also obtained using the standard quantum mechanic 
perturbation theory, applied to mass matrix eigenstates problem (see e.g. [^), but in this case it is 
very difficult to get the simple closed expressions. 
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ica symbolic manipulation language and designed to perform automatically the MI 
expansion of QFT amplitudes evaluated in mass eigenstates basis, provided that they 
are coded using a specihc format which could be parsed by the MassToMI routines. The 
results are given in terms of MI powers and the so-called divided differences [sj of the 
loop functions. The package is able to expand any type of amplitude for Hermitian 
(scalar or vector) or general complex (fermion) mass matrices, to any requested MI 
order. 

The paper is organized as follows. In Sec. [2] we formulate the algorithm used for the 
expansion. In Sec. [3]we present the syntax for dehning the amplitudes in the MassToMI 
package, the routines provided for the users and the output format. In Section 0] we 
illustrate the applications of the MassToMI package with several examples, hnally we 
conclude in Sec. 01 The MassToMI Mathematica code can be downloaded from the 
address 


WWW.fuw.edu.pl/masstomi 

2. Flavor and mass amplitudes in QFT calculations 

2.1. Mass matrices in physical QFT models 

The physical spectrum of any QFT model is dehned by the structure of the quadratic 
terms in the Lagrangian. Assuming that quadratic kinetic (momentum-dependent) 
terms have been transformed to the canonical form by appropriate held and coupling 
redehnitions and, if necessary, the Spontaneous Symmetry Breaking (SSB) mechanism 
has been used to identify physical degrees of freedom, three type of mass matrices can 
appear in the Lagrangian: 

1. Hermitian (squared) mass matrices for scalar and vectors helds, M| = (M|)l, 
diagonalized by a unitary transformation Z: 

Z^M|Z = m| =diag{ml,...,ml). (2.1) 

2. General complex mass matrices for Dirac fermions Md, diagonalized by two uni¬ 
tary transformations U, V: 

V^MdU = niD = diag(mi,..., m„). (2.2) 

The matrices V and U diagonalize also the Hermitian matrices MdMq and 
MqMd, through the transformations 

V+MdM^V = U^MijMDU = m^. (2.3) 

3. Symmetric complex mass matrices for Majorana fermions, Mn = In such 

case one can assume u = V* = o in eq. fl2.3p . so that the mass matrix is diago¬ 
nalized by a single unitary transformation O: 

O'^MnO = niN = diag(mi,... ,m„). (2.4) 

The matrix O diagonalizes also the Hermitian matrix 

OtMj,,MNO = m^. (2.5) 
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Note that the squared mass matrices of physical particles, i.e. M|, Md Md 

and M|,^ Mn must be (semi-) positive-dehnite for well-dehned QFT theories. 


2.2. Structure of the transition amplitudes 

Applying the rotations Z,U,V,0 to the held multiplets one gets the Lagrangian 
of the theory in the basis of physical (mass eigenstates) helds. The tree level vertices 
and Feynman rules in such basis depend on the elements of mixing matrices and on the 
physical particle masses. Consequently, all transition amplitudes (to any loop order) are 
linear combinations of products of mixing matrices and propagators, or loop integrals 
being the functions of physical particle masses. The dependence of such amplitudes 
on the initial Lagrangian parameters is very complicated and non-linear. Thus, as 
mentioned in the introduction, although they are compact in form and better suited 
for numerical computations, it is often practical to use the Mass Insertion expansion 
to recover an approximate, but simpler direct dependence on the interaction basis 
Lagrangian parameters. 

Such an expansion can be done with the use of Flavor Expansion Theorem, for¬ 
mulated and proven in Ref. [l|. As argued in [l|, the mixing matrices of the internal 
particles of Feynman diagrams can appear in amplitudes only in some specihc combi¬ 
nations, namely 

Scalars, vector bosons : 

Dirac fermions : llBi /(m-) VBi f{rnl) , 

UB^ VBi , 

Majorana fermions : OBif{mf)0%, OBimif{mf)OAi, (2.6) 

where f{mj) represents symbolically the dependence of the amplitude on the physical 
masses, at tree or loop level. In the expression above, and similarly for other equations 
in the rest of the paper, we assume a summation over the repeating indices (even if 
they appear more than twice). 

Assuming that the function / is analytical, the combinations listed above can be 
formally expressed as matrix elements being functions of the squared mass matrices: 


Scalars, vector bosons : 

Dirac fermions : 
UB^f{ml)Ul, 

VBi mif{m‘f)UXi 
UBi mif{ml)VXi 


f{Ml)BA , 

/(Mi^MD)^^ , 
/(MdMJ5)ba , 

\ / BA 

V / BA 


(2.7) 

(/(MDMt,)MD) , 

V / BA 

(/(Mt,MD)ML ,(2.8) 

V / BA 
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Majorana fermions : 

OB^f{m^)0% = f{Ml,M^)BA = f{M^Ml,)AB , 

OB^mJ{mi)OAi = (Mj^/(MNMj^)) = , 

V / BA \ /BA 

= (Mn/(mJ,Mn)) = (/(MnMJ,) Mn) • (2.9) 

V / BA \ /BA 

RHS of all of the expressions above depends on the matrix elements of functions of the 
hermitian matrices. As proven in Ref. [ij, such elements can be expanded in terms of 
“divided differences” of the function /, dehned recursively as 


/^(x) 


/f^+^'(xo,...,Xfc,Xfc+i) 


= fix) , 

^ fjxo) - fjxi) 

~ Xq-Xi 


/W(Xo,...,,Xfc-i,Xfc) - f^^Kxo,---,Xk-l,Xk+l) ^^ 


or, for degenerate case (and analytical /), 

1 

..- 

The expansion can be done by splitting the mass matrices into diagonal parts and 
off-diagonal mass insertions, 

+ M^ (2.12) 


where, by dehnition, 

(MS), = Mf, , 

M'fj = Ml„ Mf, = 0, (/, J = 1,..., n) . 

Then, for any analytical Hermitian matrix function /(M^) (see [ij for detailed 
tions), a matrix element will be given by the expansion 


(2.13) 

assump- 


fiM^)u = 6ufiiM^o)i) +f^\iM^o)i,iM^o)j)Mh 

(2.14) 

Ki 

+ Y. (Vot J. {/<)k„ {Mi)K,) MhM.K.Mkj + ■ ■ ■ . 

Ki,K2 


Notice that due to the dehnition 02.131) terms proportional to Mjj will always vanish 
in the summation. 
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2.3. Ambiguity in expansion of fermionic amplitudes 

It is important to note that although the alternative forms of the RHS in the last 2 
lines of eqs. fl2.8p and fl2.9p . appearing in expansion of fermionic amplitudes, are order- 
by-order equivalent when the function / is expanded in a Taylor series, they are not 
order-by-order equivalent when FET expansion of eq. fl2.14p for / is used. For example, 
lets consider the lowest order expansion of the two forms of the RHS in the last line 
of eq. fl2.9p . It would lead to (for Majorana fermions Mn = M^): 

or (2.15) 

0%,mj{mj)0% = (/(MnMJ„)Mn)^^ ^ (Mn)ba / 

with both forms obviously different for A ^ B. 

The correct choice of the form of the RHS of eqs. fl2.8p and fl2.9p and minimal 
order of expansion reproducing the exact result with sufficient accuracy depends on a 
given amplitude. To illustrate it with a realistic example, lets consider the chargino 
contribution to the lepton self energy in the Minimal Supersymmetric Standard Model 
(MSSM), illustrated in Fig. [T] 

vk 

/«' N 

/ \ 

Ij -- Ij 

Xi 

Figure 1: Sneutrino-chargino loop contributing to the lepton self-energy in the MSSM. /, J denote 
external lepton flavors. 


For the purpose of our example lets consider only the scalar part of the lepton 
self-energy, where, neglecting the small external momenta. 


Sli = SA,' 5^ my ^ g,Y‘ ^^(2.16) 

K=l i=l K=1 

where bo is the 2-point loop function, by Y) we denote the lepton Yukawa coupling, Zy 
is the unitary matrix diagonalizing the sneutrino mass matrix and Vc, Uc diagonalize 
the chargino mass matrix (for details of notation see Refs. [^, E^): 


M 2 ^ 

92V1 

V2 


Vci ‘bI f lCc = V<iAfcCc = diag(m,.,m,,). 


(2,17) 


In addition by Fk we denoted the factors in eq. 02.1611 depending only on chargino 
mixing matrices: 

2 

Fk = E ■ ( 2 - 18 ) 

i=l 
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Expanding Fk to the first MI order, one gets two possible expressions, corresponding 
to two forms of the RHS of eq. fl2.8jl . They read either as 





where we used the relation connecting divided difference of the 2-point loop function 
with the 3-point loop function, = Co{ml,ml,ml). 

In the lowest MI order arguments of bo function in eqs. fl2.19p and fl2.20p differ, which 
could be numerically important if ji and M 2 are splitted signihcantly comparing to some 
average scale Msusy- However, using the identity (holding for vanishing external loop 
momenta): 

bo{ml, ml) + ml co(ml, ml, ml) = bo{ml, ml) ml Coimj, ml, ml) , ( 2 . 21 ) 

one can see that in the 1st MI order eqs. 02.191) and 02.201) differ only by terms of the 
order of 9lvl(^2)/MsusY /Mg^gY (so of the higher order in MI expansion), in 

practice usually negligible in calculations of the leptonic transitions in the MSSM. 

3. Expanding amplitudes with MassToMI package 

Eqs. O 2 . 7 p - O 2 . 9 i) and 02.14p allow to expand any transition amplitude calculated in 
the mass eigenstates basis as a series in mass insertions powers. However, for realistic 
models with many helds it may lead to lengthy expressions, requiring tedious and error- 
prone calculations. To facilitate the problem, the MassToMI package automatizes such 
calculations. 

3.1. MassToMI installation 

The MassToMI package does not require any special installation procedures. It 
should be unpacked to the directory accessible to Mathematica, or installed system- 
wide using Install command from Mathematica menu. Then the package can be loaded 



























using the commands 
Needs["MassToMI'; 
or 

« MassToMI.m 

3.2. Syntax for constructing amplitudes in MassToMI 

In general, any amplitude, IPI irreducible or reducible, can be constructed as a sum 
of products of mixing matrices of scalars, fermions or vector bosons, physical particle 
masses, loop integrals and other factors which are not affected by the FET expansion. 
To identify the combinations of objects relevant for the expansion procedure, MassToMI 
package require them to be denoted using the special syntax, listed in Table [U (other 
symbolic or numerical factors in the amplitude, not listed in Table [H are treated by 
program as constants). 

The sub-expressions in the amplitude which could be expanded in terms of Mis need 
to have one of the forms listed on the LHS of eqs. fl2.7p - fl2.9p . Some further remarks 
are in order here: 

• As in the analytical expressions of eqs. fl2.7p - fl2.9p . it is assumed that indices of 
particles appearing more than once in the amplitude, written as Mathematica 
code, are summed over without the need of specifying the sum explicitly. 

• The amplitude does not need to be IPI, i.e. products of loop functions are allowed. 

• Only linear fermion mass powers (or multi-linear for amplitudes involving more 
than one fermion) can appear as factors in the amplitude, all even powers of 
masses should be included in the dehnitions of the loop functions. 

• Eq. (12.1 op dehnes the divided differences for functions of one variable. The gener¬ 
alization to many variables is obvious: one should calculate the divided differences 
of the required order with respect to each variable separately. However, compli¬ 
cations arise when some of the particles in the loop have equal masses (as for 
example in the case of flavor-diagonal photon couplings) and several arguments 
of loop function are identical. Then, the generalization of eq. 12.101 is less trivial 
and leads to some form of “Leibniz-like” rule, which combinatorial complication 
grows quickly with the number of identical arguments and the order of divided 
difference. MassToMI package does not allow for repeating arguments of the loop 
functions, assuming if necessary they have been appropriately renamed, e.g.: 

ct){rn\,rn\,rn\) ^ ct]new{rn\,rn\) (3.1) 

To give an example, lets consider the 3-point scalar-fermion-fermion 1-loop ampli¬ 
tude, the triangle diagram with Dirac fermions Majorana fermions Nj and scalars 
Di circulating in the loop. In general such an amplitude is the sum of terms depending 
on mixing matrices and 3-point loop function, one of them could look like 

Amplitude = Zj;Z^^*O^^Oif*V^^*Uc''mc„C 2 i{p, q, m^., m^.) (3.2) 
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Object and its syntax in MassToMI 

Arguments & examples 

Scalar (vector) mixing matrix Z 
SMIX[P,i,j] 

P - particle symbol 
i,j - particle indices 

Example: 

= SMIX[D.I,j] 

Zl/* = Conjugate[SMIX[U,I,j]] 

Left Dirac fermion mixing matrix V 

FMIXL[P,i,j] 

P - particle symbol 
i,j - particle indices 

Example: 

= FMIXL[D.I,j] 

Right Dirac fermion mixing matrix [/ 

FMIXR[P,i,j] 

P - particle symbol 
- particle indices 

Example: 

Uq* = Conjugate [FMIXREQ, I ] 

Majorana fermion mixing matrix 0 
NMIX[P,i,j] 

P - particle symbol 
i,j - particle indices 

Example: 

Oj' = NMIX[Q,I,j] 

physical particle mass 

MASS[P,i] 

P - particle symbol 
i - particle index 

Example: 

= MASS[Q,I] 

general loop fnnction 

LOOP[name,{{Pl,il}{a,...}] 

name - name of function 
{{PI, if}, • • •} - physical loop masses 
given as the list of (particle,index) pairs 
{a,...} - optional other arguments 

Examples: 

LOOP[b21,{{U,i},{D,j}},{p}] 


L00P[c0,{{A,i},{B,j},{C,k}},{p,q}] 


Table 1: Syntax of objects used to constructing amplitudes in MassToMI package 
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To expand this in terms of mass insertions with the use of MassToMI routines, one must 
rewrite it into a Mathematica expression as: 

Amplitude = SMIX[D,I,i] Conjugate [SMIX[D,J,i]] 

NMIX [N,K,j] Conj ugate [NMIX [N,L,j]] 

Conjugate[FMIXL[C,M,n]] FMIXR[C,N,n] MASS[C,n] 
LDDP[c21,{{C,n},{D,i},{N,j}},{p,q}] ; 


3.3. Control variables in MassToMI 

MassToMI package uses the following control variables, dehning which combinations 
of mixing matrices should be expanded to which order in terms of mass insertions, and 
in which form the hnal result should be presented. 

• FetScalarList={{Pl, ol}, . . . ,{Pn,on}}. In this list the user can specify sym¬ 
bols for the scalar and vector particles (PI, . . . ,Pn) for which the mixing matrices 
are expanded in MI powers to orders ol, . . . ,on, respectively. 

• FetFermionList={{Pl , ol ,M1} , . . . , {Pn, on,Mn}}. Specihes list of Dirac or Ma- 
jorana fermions PI, . . . ,Pn for which the mixing matrices are expanded to MI 
orders ol, . . . , on. Each of the arguments Ml, . . . ,Mn can take one of two values, 
MHM or MMH, deciding which form of the RHS in the last 2 lines of eqs. fl2.8p and 
fl2.9p is used, depending on M^M or on MM^, respectively (see discussion in 
Section 12.3p . 

• FetMaxOrder=n. Only mass insertion products of the total order FetMaxOrder 
or lower are kept in the hnal result. 

Particles with symbols not specihed on FetScalarList or FetFermionList lists are 
left unexpanded. 

3 . 4 .. The main expansion routine and output syntax 

The actual MI expansion is done by a call to the function 

FetExpand[ Amplitude ] 

The FetExpand routine performs the following actions: 

1. It checks hrst the structure of the amplitude for possible syntax problems, dis¬ 
playing if necessary the relevant error or warning messages. 

2. It hnds pairs of mixing matrices for which the second index is identical (like 
e.g. in Xq^*Xq^ product), and the particle symbol {Q in this case) appears on 
FetScalarList or on FetFermionList (as mentioned above it is assumed that 
the repeating index denotes summation over it). For pairs of fermionic matrices, 
the program also searches if they are associated with any linear fermion mass 
factors with the same summation index. 
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Object 

Syntax and arguments 

Example 

Diagonal entry of scalar 
squared mass matrix 

MS2[P,i,i] 

P - particle symbol 
i - particle index 

= MS2[P.i,i] 

Fermion mass matrix 

MF [P,i,j] 

P - particle symbol 
- particle indices 

{Mp),j = MF[P,i J] 

Diagonal entry of squared M^M 
fermion mass matrix 

MHM[P,i,i] 

P - particle symbol 
i - particle index 

{Ml,Mp)ii = MHM[P,i,i] 

Diagonal entry of squared MM^ 
fermion mass matrix 

MMH[P,i,i] 

P - particle symbol 
i - particle index 

{MpMl)ii = MMH[P,i,i] 

Scalar MI matrix 

MS2I[P,i,j] 

P - particle symbol 
qj - particle indices 

= MS2I[P,iJ] 

Fermion squared M MI matrix 

MHMI[P,i,j] 

P - particle symbol 
i,j - particle index 

{Ml,Mp)ij = MHMI [P.i.j] 

Fermion squared MM^ MI matrix 

MMHI[P,i,i] 

P - particle symbol 
i,j - particle index 

{MpMl)ij = MMHI[P,i,j] 


Table 2: Notation for the mass matrices in the interaction basis appearing in the expanded amplitudes. 
The objects denoting off-diagonal mass insertions, MS2I, MHMI and MMHI, are assumed to have O’s on 
diagonal. 


3. The routine finds the loop functions which arguments have indices identical to 
those in the matching pairs of mixing matrices and performs the FET expansion 
for each of such pairs, using the formula of eqs. fl2.7jl - fl2.9jl and eq. fl2.14jl . Some 
further syntax checks are performed at this stage and reported if necessary. 

4. Finally the result is displayed in terms of the objects used to denote diagonal and 
off-diagonal entries of the mass matrices listed in Table [2l 

The expanded amplitude depends on the divided differences of the loop functions 
appearing in the original unexpanded expression. The notation for the divided differ¬ 
ences is as follows: 

• The loop functions are renamed using the first argument of the original expression; 
the new first argument (non-negative integer) is added to denote the order of 
divided difference. The list containing the other arguments is flattened, so that 
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the arguments are no longer splitted into masses and optional variables. This is 
illustrated by the following example: 

LDDP[b21,{{U,i},{Cj}},{p}] ^ b21[0,{U,i},{Cj},p ] 

• The physical mass arguments of the loop functions undergoing the MI expansion 
are replaced by a list containing the diagonal entries of mass matrices (boson or 
fermion squared masses) representing the variables of divided difference in the 
given argument. 

• The divided differences of the loop functions are multiplied by the relevant mass 
insertion factors. Again, it is assumed that all repeated indices named fetQxx 
should be summed over. 

Continuing the example given above, lets assume that MI expansion has been 
performed to 2nd order for a fermion C. Then the highest order term in the 
expansion is proportional to the function denoted as: 

Coiijugate[FMIXL[C,c,j]] FMIXL[C,d,j] LOOP[b21,{{U,i},{C,j}},{p}] 

MMHI[C,d,fetCl] MMHI[C,fetCl,c] x 

X b21[2,{U,i},{MMH[C,d,d],MMH[C,fetCl,fetCl],MMH[C,c,c]},p] 

where 2 being used as the 1st argument of the b21 function is the order of the 
divided difference, MMH [C , i, i] = ,,, MMHI [C, i»j] = (Afc-M^)for i ^ j 

and the summation over the internal index fetCl is assumed. 

• If the FET expansion is performed in several indices, the hrst argument of the 
expanded loop function denotes the total order (the sum of all orders) of the 
divided differences in all arguments. 

• If the hnal result still depends on some of the physical masses (not expanded into 
MI series), they are now denoted as MASS [P, i], e.g. 

b21[2,{U,i},{MMH[C,d,d],MMH[C,fetCl,fetCl],MMH[C,c,c]},p] 

b21 [2 ,MASS [U, i] 2,{MMH[C,d,d] ,MMH[C,fetCl,fetCl],MMH[C,c,c]},p] 

• If necessary, the divided differences can be re-expressed as the explicit combina¬ 
tions of the original loop functions with the use of FetExpandDividedDif f erences 
functions described in next section. 

3.5. Auxiliary functions 

In addition to the main FetExpand routine, the MassToMI package provides sev¬ 
eral auxiliary functions. For consistency and to distinguish them from other routines 
provided by Mathematica, all their names, by convention, start from the prehx Fet. 

One of the auxiliary functions, FetExpandDividedDif f erences, allows to re-express 
divided differences appearing in the FetExpand output as the combinations of the ini¬ 
tial loop functions. Other functions help to manipulate expressions with repeating 
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Function Arguments 


FetExpandDividedDifferences[ DivFunction, (ReplacementRule) ] 

DivFunction divided difference 

ReplacementRule rule how to redefine Fnnction (optional) 

FetSumWithDelta[ Expression, (ExcludeList) ] 

Expression expression containing Kronecker 5-symbols 

ExcludeList index or list of indices excluded from summation (optional) 


FetSumParticle[ Expression, Particle, Range, (ExcludeList) ] 

Expression expression to expand repeating indices 

Particle symbol of particle (or list of particles) which indices are expanded 

Range={Low,Up} summation range for repeating indices 

ExcludeList index or list of indices excluded from summation (optional) 


FetSumFactor[ Expression, 
Expression 
Factor 

Range={Low,Up} 
ExcludeList 


Factor, Range, (ExcludeList) ] 
expression to expand repeating indices 
factor name (or list of factors) which indices are expanded 
summation range for repeating indices 
index or list of indices excluded from summation (optional) 


FetSumIndexE Expression, Index, Range ] 

Expression expression to expand 

Index name of index to expand 

Range={Low,Up} summation range for Index 


Table 3: List and arguments of auxiliary functions provided by MassToMI package. 

indices, assumed to be implicitly summed over. They naturally appear in the higher 
order terms of the FET expansion (see eq. fl2.14p ) in sums over mass matrix indices. 
In many cases, the summation convention can be also used to dehne the initial (unex¬ 
panded) amplitude in a compact form. However, for the analysis of the physical effects 
it may be useful to expand some (or all) repeated indices as explicit sums. For that 
purpose, MassToMI package provides 3 routines, FetSumParticle, FetSumFactor and 
FetSumIndex. In addition, Mathematica, at least in version 10, does not have summa¬ 
tion convention rules with Kronecker (5-symbol implemented. Resumming Kronecker 
symbols can be done using the FetSumWithDelta function. 

Below we describe in details the syntax for using the auxiliary functions (their 
compact list is collected in Table |3]). One should note that some of their arguments, 
given in parenthesis with slanted font, are optional. 

3.5.1. FetExpandDividedDifferences[ DivFunction, (ReplacementRule) ] 
FetExpandDividedDif f erences expands the divided difference of a function given as 
DivFunction = fname[(n + ..., {M_l ,...,M_(n+l)},...,optional arguments] 
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into an explicit combination of the initial fnnctions with different mass argnments (the 
hrst argnment, order of divided difference, is trnncated in the expanded form). For the 
degenerate mass arguments relevant derivatives (in the standard Mathematica nota¬ 
tion) of the initial function are used. Examples are given below. 

Non-degenerated arguments: 

FetExpandDividedDifferences[ b21[2,{ml,m2,m3},M1,p] ] = 

b21[ml,Ml,p]/(ml - m2)/(ml - m3) - b21[m2,Ml,p]/(ml - m2)/(m2 - m3) 

+ b21[m3,Ml,p]/(ml - m3)/(m2 - m3) 

Degenerated arguments: 

FetExpandDividedDifferences[ b21[2,{ml,m2,ml},Ml,p] ] = 

(b21[m2,Ml,p] - b21[ml,Ml,p])/(ml - m2)^ + b2l(^’°’°) [m2,M1 ,p] /(ml - m2) 

FetExpandDividedDifferences [ b21 [2,{ml ,ml ,ml} ,M1 ,p] ] = b2l(^’*^’°) [ml ,M1 ,p]/2 

Divided difference of multiple arguments, of total order 2 -|- 1 = 3: 
FetExpandDividedDifferences [ b21[3,{ml,ml,ml},{M1,M2},p] ] = 

(b2i(2,0,0) [nii^Ml.p] - b2l(2A0) [nii^M2,p])/2/(Ml - M2) 

The second optional argument of FetExpandDividedDifferences should be speci- 
hed if the function used as the hrst argument was renamed to avoid multiple identical 
mass arguments, as in the example given in eq. fl3.ip . Specifying the replacement rule 
as the second argument allows to revert the redehnition and to go back to the original 
function name, argument list and, for the degenerate mass arguments, to calculate cor¬ 
rectly the function derivatives. Example: 

FetExpandDividedDifferences [ cOnew[2,{ml,m2,ml},M1] ] = 

(cOnew[m2,M1] - cOnew[ml,M1])/(ml - m2)^ 

+ cOnewd’*^) [ml ,M1] / (ml - m2) 

but 

FetExpandDividedDifferences [ cOnew[2,{ml,m2,ml},M1], cOnew[a_,b_]->c0[a,b,a] ] = 
(cO[m2,Ml,m2] - cO[ml,M1,ml])/(ml - m2)^ 

+ (cOd-0-0) [ml,Ml,ml] + cO(°’°,i) [mi, m 1 ,ml] )/(ml - m2) 

3.5.2. FetSumWithDelta[ Expression, {ExcludeList) ] 

FetSumWithDelta function searches for occurrences of the factor KroneckerDelta[a,b] 
in Expression and if possible replaces them using the standard rule 

KroneckerDelta[a,b] X[a] -> X[b], 

excluding the cases when summation index is specihed on the (optional) ExcludeList. 
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This option allows e.g. to exclude external particle indices from summation. Example: 
exp = KroneckerDelta[I,a] KroneckerDelta[J,b] Fun[a,b] 

FetSumWithDeltaf exp ] = Fun[I,J] 

FetSumWithDeltaf exp, {b} ] = KroneckerDelta[J,b] Fun[I,b] 

Few additional remarks are in order: 

• FetSumWithDelta does not perform summation if the repeating index is not a 
simple variable but an expression by itself, so terms of the form 

KroneckerDelta[a+1,b] f[a] 

are left unchanged. However the argument of f is evaluated properly: 
FetSumWithDelta[ KroneckerDelta[a,b] f[a+1] ] = f[b+1]. 

• If the second KroneckerDelta argument, complementary to the summation index, 
is an expression, replacement is done: 

FetSumWithDelta[ KroneckerDelta[a,b+3] f[a] ] = f [b+3] 

However, FetSumWithDelta does not check for the allowed range for summation 
indices - if in the example above a,b = 1. . .3, then KroneckerDelta [a, b+3] = 
0 and the result given by FetSumWithDelta is not correct. Expressions of that 
type must be simplihed manually by the user. 

3.5.3. FetSumParticle[ Expression, Particle, Range, (ExcludeList) ] 

FetSumParticle expands the summation convention for the repeating indices of 
mass insertion matrices for given set of particles specihed in the argument Particle (it 
could be a single particle symbol or list of symbols). More specihcally, it searches for 
the occurrence of factors of the form X [Particle, a,b], where X is one of the following 
matrices: X=MF,MS21 ,MHMIIf such a factor is found, and at least one of its 
indices (not specihed on the optional argument ExcludeList) is present also in the 
term multiplying it, the repeating index is replaced by the explicit sum in the range 
Range={Low,Up}. Examples: 

exp = MF[Q,B,K] MF[Q,B,J] + MS2I[P,K,A] MA[A,J] 

FetSumParticle [ exp,P,{1,2} ] = ]yiF[Q,B,J] MF[Q,B,K] 

+ MA[1,J] MS2I[P,K,1] + MA[2,J] MS2I [P,K,2] 

FetSumParticle [ exp,{P,Q},{l,2} ] = MF[Q,1,J] MF[Q,1,K] 

+ MF[Q,2,J] MF[Q,2,K] + MA[1,J] MS2I [P, K, 1] + MA[2,J] MS2I[P,K,2] 

Note that whenever possible, FetSumParticle (and also other functions dealing with 
summation convention described below) automatically sets to zero diagonal elements 
of the mass insertion matrices, like MS21 [P,2,2] factor in the example below: 

FetSumParticle [ exp /. K->2 ,{P,Q},{1,2} ] = MFfQ,!,!] MF[Q,1,2] 

+ MF[Q,2,J] MF[Q,2,2] + MA[1,J] MS2I [P,2,1] 
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3.5.4- FetSumFactor[ Expression, Factor, Range, (ExcludeList) ] 

FetSumFactor routine can be used to expand summation over indices of objects not 
related to MI expansion (like CKM matrix). It searches for the repeating indices of fac¬ 
tor (s) specihed as the argument Factor. Factors can have any number of indices, but 
they are all summed in the same Range={Low,Up}. If various indices have different sum¬ 
mation ranges, routine should be called several times, specifying optional ExcludeList 
argument to avoid summation over indices not belonging to correct Range. Examples: 

exp = MA[A,K] CKM[A,J] 

FetSumFactor[ exp,CKM,{l,3} ] = CKM[1,J] MA[1,K] + CKM[2,J] MA[2,K] 

+ CKM[3,J] MA[3,K] 

FetSumFactor[ exp,CKM,{l,3}, A ] = CKM[A,J] MA[A,K] 

3.5.5. FetSumIndex[ Expression, Index, Range ] 

FetSumIndex routine can be used to expand summation over single index directly 
specihed by user. Example: 

exp = 1 + A[J,K] B[J,L] 

FetSumIndexE exp,J,{1,3} ] = 1 + A[1,K] B[1,L] + A[2,K] B[2,L] + A[3,K] B[3,L] 


3.5.6. Special cases and limitations 

Syntax of MassToMI package assumes in general that repeating indices are implicitly 
summed over. However, this is not a summation convention in a “classical” sense, 
as it allows for indices repeating more than twice. Such dehnition is convenient for 
applications of FET expansion, but, used without proper care, may lead to various 
ambiguities. Thus, it is advisable to check the structure of the intermediate expressions 
appearing during calculations, particularly before applying to them auxiliary functions 
designed for manipulating objects with repeating indices. 

Few additional remarks may be helpful to avoid problems with the incorrect usage 
of routines provided by MassToMI package: 

• It is highly recommended to use the ExcludeList argument of FetSumWithDelta, 
FetSumParticle and FetSumFactor functions for specifying which indices cor¬ 
respond to external particles in given amplitude, to avoid resumming over them 
accidentally. 

• Notice that the functions FetSumFactor, FetSumParticle and FetSumWithDelta 
search only for repeating indices in simple factors in the analyzed expression, not 
checking for their appearance in lower level sub-expression (those are left un¬ 
changed). For example: 

FetSumWithDelta[ KroneckerDelta[a,b] f[a] ] = f[b] 
but 
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FetSumWithDeltaE Fun[KroneckerDelta[a,b] f[a]] ] = 

Fun[f[a] KroneckerDelta[a,b]] 

and similarly for FetSumParticle and FetSumFactor. FetSumIndex will work 
even for indices hidden in sub-expressions, assuming that they appear in at least 
two different factors. 

• Functions FetSumParticle and FetSumFactor treat higher powers of matrix ele¬ 
ments as products with repeating indices, and perform appropriate resummations, 
as illustrated below: 

FetSumFactor [ M[a,b]''2,M,{1,2} ] = FetSumFactor [ ]yi[a,b]]yi[a,b] ,]yi,{l,2} ] = 

+ M[1,2]2 + M[2,1]2 + M[2,2]2 

• One should strictly avoid using the same names for various types of quantities 
appearing in the transition amplitudes or other expressions used as arguments of 
MassToMI functions. In particular, functions dealing with summation convention 
cannot distinguish what is really an “index” and try to sum over all symbols in 
given expression with the names identical to repeating indices. This may lead to 
strange bugs and obviously incorrect or nonsensical results, like in the examples 
given below: 

FetSumWithDeltaE KroneckerDelta[A,B] A[A] ] = B[B] 

FetSumFactorE M[A,B] A[A],A,{1,2} ] = 1[1] M[1,B] + 2[2] M[2,B] 


4. Examples of the MassToMI applications 

We collect below several examples illustrating the functionality of the MassToMI 
package, from simple test amplitudes to the realistic case of the expansion of one of 
the diagrams contributing to the flavor violating decay of the Higgs boson to leptons, 
h r/i, in the MSSM. 

4 . 1 . Expansion of simple amplitudes 

The two basic examples of MassToMI usage presented below are included in the file 
mmi_example .m attached to the MassToMI distribution. 

Example 1. The simplest case involves a pair of mixing matrices of a single scalar 
particle multiplying loop function a{m?). Lets assume that amplitude should be ex¬ 
panded to 2nd order in MI powers. 

Ampl = Conjugate [SMIXfP,a,i]] SMIX[P,b,i] LDDP[A0,{{P,i}}]; 

FetScalarList = {{P,2}}; 

FetMaxOrder = 2; 

FetAmpl = FetExpand[Ampl]; 

Print["Expanded amplitude = ", FetAmpl]; 
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The result is: 

Expanded amplitude = AO[0,MS2[P,b,b]] KroneckerDelta[a,b] + 

A0[l,{MS2[P,b,b] ,MS2[P,a,a]}] MS2I[P,b,a] + 

A0[2,{MS2[P,b,b],MS2[P,fetPl,fetPl],MS2[P,a,a]}] MS2I [P,b,fetPl] MS2I[P.fetPl,a] 

Note that only amplitudes of a specihc structure (see eqs. fl2.7l) - fl2.9p l are allowed and 
FetExpand reports errors if incorrect combinations of mixing matrices are used. For 
example, if c.c. is removed from one of the scalar matrices in the example above, 

Ampl = SMIX[P,a,i] SMIXLP.b.i] LDDP[A0,{{P,i}}]; 
program reports: 

MassToMI::cmplx2: ERROR: Incorrect matching of c.c.’s in product of SMIX[P,a,i] 
and SMIX[P,b,i]? 

and aborts the execution. 

In other cases only warnings or informational messages are displayed but program 
still makes an attempt to calculate the result, like for instance the mismatch in the 
summation index of scalar matrices: 

Ampl = Conjugate[SMIXfP,a,j]] SMIX[P,b,i] L00P[A0,{{P,i}}]; 

leads to 

MassToMI::war1: WARNING: unmatched index of scalar on FetExpandList in 
SMIX[P,a, j]? 

MassToMI::war1: WARNING: unmatched index of scalar on FetExpandList in 
SMIX[P,b,i]? 

Expanded amplitude = AO[0,MASS [P,i]Conjugate[SMIX[P,a,j]] SMIX[P,b,i] 

meaning that no pair of scalar mixing matrices with matching 2nd index was found, so 
the result is simply equal to the initial amplitude rewritten using the syntax rules for 
the MassToMI output. 

The full list of errors, warnings and informational messages can be found in the 
header of the MassToMI .m hie. 


Example 2. Let us now consider a more complicated case of an amplitude depending 
on scalar P and Dirac fermion F masses. We request 1st order of MI expansion for 
both particles and set FetMaxOrder = 1 as a maximum total MI order. The preferred 
form of the squared fermion mass in output is M^M. 

Ampl = Conjugate[SMIX[P,a,i]] SMIX[P,b,i] Conjugate[FMIXL[F,c,j]] FMIXR[F,d,j] 
MASS[F,j] L00P[B0,{{P,i},{F,j}}] ; 

FetScalarList = {{P,l}}; 

FetFermionList = {{F,1,MHM}}; 
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FetMaxDrder = 1; 

FetAmpl = FetExpand[Ampl ]; 

Print["Expanded amplitude = ", FetAmpl]; 

The result is 

Expanded amplitude = BO[0,MS2[P,b,b], MHM[F,d,d]] Conjugate[MF[F,c,fetFl]] 
KroneckerDelta[a,b] KroneckerDelta[d,fetFl] 

+ B0[l,]yiS2[P,b,b] ,{MHM[F,d,d] ,MHM[F,fetFl,fetFl] }] Conjugate [MF[F,c,fetFl]] 
KroneckerDelta[a,b] MHMI[F,d,fetFl] 

+ BO[1,{MS2[P,b,b],MS2[P,a,a]},MHM[F,d,d]] Conjugate [MF [F,c,fetFl]] 
KroneckerDelta[d,fetFl] MS2I[P,b,a] 

where again summation over the repeating index f etFl is assumed. The result contains 
some Kronecker 5-symbols with fermion indices which can be resummed 

FetAmpl = FetSumWithDelta[FetAmpl,{a,b,c,d}] ; 
giving the simpler expression 

Expanded amplitude = BO[0,MS2[P,b,b],MHM[F,d,d]] Conjugate[MF[F,c,d]] 
KroneckerDelta[a,b] 

+ B0[l,MS2[P,b,b],{MHM[F,d,d],MHM[F,fetFl,fetFl]}] Conjugate[MF[F,c,fetFl]] 
KroneckerDelta[a,b] MHMI[F,d,fetFl] 

+ B0[l,{MS2[P,b,b],MS2[P,a,a]},MHM[F,d,d]] Conjugate[MF[F,c,d]] MS2I[P,b,a] 


4-3. h ^ Tfi decay in the MSSM. 

As a more realistic and complicated example we consider the flavor violating decay 
of the Higgs teson to leptons, h —)■ r/r, in the MSSM. This decay was recently measured 


by ATLAS llj] and CMS [l2| Collaborations, and reported to be more frequent than 
predicted within the Standard Model. It was also theoretically recently reanalyzed 
within the MSSM [l3| and claimed to be enhanced. 

For simplicity, for the purpose of our example we expand the contribution from one 
diagram only, shown in Fig. [2] (results for MI expansion of sum of all diagrams contribut¬ 
ing to the effective Higgs-fermion vertex in the MSSM can be found in Refs. (l4. 15. 16|]. 
It depends on the scalar (slepton) and fermionic (neutralino) mixing matrices. 

In the approximation of the vanishing external momenta, the amplitude correspond¬ 
ing to this diagram reads as: 

(4.1) 


A = 


I -pABK p 


R ■ 


where 


t^ABK 

^L{R) 


- C'o, (4.2) 

i,j,l 


and A, B are lepton generation indices {B = 3 and A = 2 for r/i in the final state, K 
is the Higgs index, K = 2 for the “little h” Higgs boson and iC = 1 for the “big H” 
and Cq is the standard 3-point scalar loop integral for vanishing external momenta. 
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Figure 2: Slepton-neutralino diagram contributing to the h ^ decay in the MSSM 


Following strictly the notation and conventions of Refs, j^, [l^ (we do not repeat 
them here explicitly), the Higgs-slepton and lepton-slepton-nentralino vertices read as 
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(4.3) 

(4.4) 


(4.5) 


Neglecting the small terms proportional to lepton Yukawa couplings and, for simplicity, 
the “non-holomorphic” + trilinear soft terms, the vertices and amplitude for are 

coded as (summation over repeating generation indices C, D in HLL is assumedjj; 

HLL[K_,l_,i_] = - e‘2/2/cw~2 (vl ZR[1,K] - v2 ZR[2,K]) (KroneckerDelta[i,1] 

+ (1 - 4 sw'' 2 )/ 2 /sw ''2 Conjugate [SMIX[L,C,i]] SMIX[L,C,1]) 

- 1/Sqrt[2] ZR[1,K] (Conjugate[AL[C,D]] SMIX[L,C,1] Conjugate[SMIX[L,D+3,i]] + 
AL[C,D] Conjugate[SMIX[L,C,i]] SMIX[L,D+3,1]); 

LSNL[J_,j_,i_] = e/(Sqrt[2] sw cw) SMIX[L,J,j] (NMIX[N,l,i] sw + NMIX[N,2,i] cw) ; 
LSNR[J_, j_,i_] =- e Sqrt [2]/cw SMIX[L, J+3, j] Conjugate [NMIX[N, 1,i]] ; 


^The Mathematica code for the example discussed in this section is attached to the MassToMI 
distribution as the file htaumu_decay .m 
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(* amplitude *) 

ampl = - MASS[N,j] HLL[K,l,i] Conjugate[LSNR[A,1,j]] LSNL[B,i,j] * 
LDDP[cO,{{L,i},{L,l},{N,j} }] // Expand; 

(* simplification: sum over KroneckerDelta[i,1] *) 
ampl = FetSumWithDelta[ampl,{A,B,K}] ; 

(* give new name for loop function with repeating arguments *) 

ampl = ampl /. LOOP[cO,{{L,i_},{L,i_},{N, j}}] -> LOOP [cOsq, {{L, i} , {N, j}}] ; 

At this stage the amplitude is correctly dehned in terms of MassToMI objects (expres¬ 
sions for can be obtained by replacing L ■H- R). For our purpose we expand it 

up to the lowest non-trivial order, i.e. 1st order in the slepton mass insertion and 0th 
order in the neutralino mass insertion (neutralino mass matrix entries are not related 
to flavor violation, so neglecting higher orders is equivalent to skipping 
suppressed terms). The actual expansion is done by: 

(* control variables *) 

FetMaxOrder = 1; 

FetScalarList = {{L,l}}; 

FetFermionList = {{N,0,MHM}}; 

(* MAIN EXPANSION ROUTINE *) 
fetampl = FetExpandf ampl ]; 

The expanded amplitude is a lengthy expression build from objects dehned in Sec¬ 
tion 13.41 It could be further signihcantly simplihed by substituting the explicit form 
for the neutralino and slepton mass matrices. 

(* simplify KroneckerDelta terms, excluding external indices *) 

fetampl = fetampl /. KroneckerDelta [3 + a_, 3 + b_] -> KroneckerDelta [a, b] ; 

fetampl = fetampl /. KroneckerDelta[3 + a_, b_] -> 0; 

fetampl = fetampl /. KroneckerDelta[a_, 3 + b_] -> 0; 

fetampl = FetSumWithDelta[ fetampl, {A,B,K} ]; 

(* perform explicit sum over AL indices (expand summation convention) *) 
fetampl = FetSumFactor[ fetampl, AL, {l, 3} ]; 

(* Neutralino mass matrix MN and MHMN = MN~+MN *) 

MN = Table[0,{i,l,4},{j,l,4}]; 

MN[[1,1]] = Ml; 

MN[[2,2]] = M2; 

MN[[3,4]] = MN[[4,3]] = mu; 

MN[[1,3]] = MN[[3,1]] = - e vl/2/cw; 

MN[[1,4]] = MN[[4,1]] = e v2/2/cw; 

MN[[2,3]] = MN[[3,2]] = e vl/2/sw; 

MN[[2,4]] = MN[[4,2]] = - e v2/2/sw; 
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MHMN = ConjugateTranspose[MN].MN; 


(* Slepton mass matrix ML2 (neglecting small Yukawa contributions) *) 

ML2 = Table[0,{i,l,6},{j,l,6}]; 

For [i=l,i<4,i++, For[j =1,j <4,j ++, 

ML2[[i,j]] = e''2 (vl~2-v2~2) (1-2 cw~2)/(8 sw''2 cw~2) KroneckerDelta[i,j] + MLL[j,i] 
ML2[[i+3,j+3]] = - e~2 (vl''2-v2''2) /(4 cw"2) KroneckerDelta[i,j] + MRR[i,j]; 
ML2[[i,j+3]] = vl/Sqrt[2] AL[i,j]; 

ML2[[i+3,j]] = Conjugate[ML2 [ [i,j+3]]]; 

] ]; 

(* slepton MI matrix *) 

ML2I = ML2; 

For[i=l,i<7,i++, ML2I[[i,i]]=0 ]; 

(* choose indices for h->tau mu decay *) 
fetampl = fetampl /. A->2 /. B->3 /. K->2; 

(* substitute real form of neutralino and slepton mass matrices *) 

For [i=l,i<5,i++, For[j =1,j <5,j ++, 

fetampl = fetampl /. MF[N,i,j] -> MN[[i,j]] /. MHM[N,i,j] -> MHMN[[i,j]]; 

]; ]; 

For [i=l,i<7,i++, For[j =1,j <7,j ++, 

fetampl = fetampl /. MS2[L,i,j] -> ML2[[i,j]] /. MS2I [L,i,j] -> ML2I[[i,j]]; 

]; ]; 

(* remove unnecessary c.c. from real parameters *) 
fetampl = fetampl // FunctionExpand; 

fetampl = fetampl /. Conjugate [e] -> e /. Conjugate [cw] -> cw 

/. Conjugate[sw] -> sw /. Conjugate[vl] -> vl 
/. Conjugate[v2] -> v2 // Expand; 

(* for simplicity assume real Ml *) 
fetampl = fetampl /. Conjugate[Ml] -> Ml; 

At this stage the flavor violating terms could come from the FET expansion or directly 
from the trilinear Ai terms explicitly present in the Higgs-slepton vertex of eq. fl4.5p , so 
their product could be of a order higher then 1. In addition, leaving the (9(M^/M|^gy) 
terms in the arguments of loop functions is inconsistent, as similar terms have been 
neglected in expanding neutralino mass matrices to 0th order only. Thus, further sim¬ 
plifications should be done: 

(* kill higher order MI terms *) 

fetampl = fetampl /. AL[A_,B_] -> eps AL[A,B] /. ALP[A_,B_] -> eps ALP[A,B] /. 
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MLL[A_,B_] -> eps MLL[A,B] /. MRR[A_,B_] -> eps MRR[A,B]; 

fetampl = fetampl /. AL[A_,A_] -> AL[A,A]/eps /. ALP[A_,A_] -> ALP [A, A]/eps / 

MLL[A_,A_] -> MLL[A,A]/eps /. MRR[A_,A_] -> MRR[A,A]/eps; 
fetampl = Normal[Series[fetampl,{eps,0,l}]] /. eps->l; 

(* neglect terms 0(]yiW~2/]yiSUSY~2) in loop function arguments *) 
fetampl = Simplify [fetampl] /. e~2(vl~2 + v2''2) -> 0 

/. e''2(vl~2 - v2*2) -> 0 /. e*2(-vl''2 + v2*2) -> 0; 

(* define and substitute explicit Higgs mixing matrix and vev’s *) 

ZH = Table[0,{i,l,2},{j,l,2}]; 

ZH[[1,1]] = ZH[[2,2]] = Cos [alpha]; 

ZH[[1,2]] = - Sin[alpha]; 

ZH[[2,1]] = Sin[alpha]; 
vl = 2 MW sw/e Cos[beta]; 
v2 = 2 MW sw/e Sin[beta]; 

For[i=l,i<3,i++. For[j=l,j<3,j++, fetampl = fetampl /. ZR[i,j]->ZH[[i,j]] ] ] 

In spite of lengthy and complicated intermediate expressions, difficnlt to obtain without 
the use of computer, the hnal result is compact and simple: 

fetampl = 3 e^ Ml/(2 Sqrt [2] cw^) ( 

+ 2/3 cw2 AL[2,2] cO [1 ,MLL [3,3] ,MLL [2,2] ,MRR [2,2] .Ml^] MLL[2,3] Sin[alpha] 

+ 2/3 cw2 AL[3,3] cO [1 ,MLL [3,3] ,MRR [3,3] ,MRR[2,2] .Ml^] MRR[3,2] Sin[alpha] 

+ (2 cw^ c0[0,MLL[3,3] ,MRR[2,2] ,M12] Sin[alpha] 

+ 2 MW^ ((4 cw^ - 3) cO[l,MLL[3,3] ,MLL[3,3] ,MRR[2,2] ,M12] 

+ 2 sw2 c0sq[l,MLL[3,3],MRR[2,2],M12]) Cos [beta] Sin[alpha+beta])/3 AL[3,2] ) 

As expected, it is linear in the three flavor violating slepton mass insertions, MLL [2,3] , 
MRR [2,3] and AL [3,2], with the coefficients given in a simple analytical form. 

The divided differences of loop functions can be further expanded with the use of 
FetExpandDividedDif f erences routine. The latter step may be in particular required 
if the diagonal slepton mass terms are degenerated and divided differences should be 
replaced by relevant derivatives. This may require some care: e.g. the cOsq function 
has been earlier dehned in the code as c0sq[ml^, m2^] = c0[ml^, ml^, m2^], so the divided 
differences expansion should be done reverting the redehnition: 

FetExpandDividedDifferences [c0sq[ 1, {MLL[3,3],MRR[2,2]},M1^] , 

c0sq[a_,b_] -> c0[a,a,b] ] = 

(cO[MLL[3,3] ,MLL[3,3] ,M12] - cO [MRR [2,2] ,MRR [2,2] ,Ml2] ) / (MLL [3,3] -MRR[2,2]) 

Obviously, the same program could be easily adapted to expand other diagrams for this 
process and hnally to calculate the decay branching ratio directly in terms of parameters 
of the initial MSSM Lagrangian. If necessary, with minimal modihcations, it could be 
also used to obtain higher order terms, both in flavor violating mass insertions and in 
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M'^/Mgjjgy powers. Such higher order terms may appear important, especially when 
lower order MI powers cancel out, as it often happens for rare decays in the MSSM - 
see e.g. discussion of t —)■ ch, uh decays in Ref. [iTj , where the accuracy of the results 
expanded by FET technique has been compared vs. the exact numerical calculations 


performed in mass eigenstates basis with the use of SUSY_FLAVOR library |18l. Il9l. |20 


5. Summary 

We presented MassToMI vl.O, a Mathematica package for an automatic expansion 
of transition amplitudes calculated in the mass eigenstates basis in terms of series of 
mass insertions. The expressions for the mass eigenstates amplitudes are usually more 
easier to calculate diagrammatically and better suited for the numerical computations 
but depend on the initial (interaction basis) Lagrangian parameters in a complicated 
way. The MassToMI routines allow to obtain an analytical approximation for the tran¬ 
sition amplitudes directly as a power series in terms of the off-diagonal entries of mass 
matrices, without the need of a separate calculation of the Feynman diagrams with 
mass insertions as vertices. 

The package is general enough to expand any amplitude, of any loop order, in any 
model involving scalar, fermion (Dirac or Majorana) or vector particles, expressing the 
result in terms of the divided differences of the loop functions. It can perform the Mass 
Insertion expansion to high orders, limited only by the combinatorial complication of 
the result and growing computational time. 

In addition, MassToMI provides several auxiliary functions allowing to express di¬ 
vided differences of any order in terms of the initial function used to generate them, 
and to manipulate and expand expressions containing factors with repeating indices, 
assumed to be implicitly summed over. 

The current version of MassToMI Mathematica code and its mannal can be down¬ 
loaded from the address 


WWW.fuw.edu.pl/masstomi 
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